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I. INTRODUCTION 



A well-known programming problem is the traveling- salesman 
problem . As classically stated, a salesman must start from his home- 
town, travel to each of n-1 other cities, and return home. He is re- 
quired to enter each of the n cities exactly once and leave each city 

exactly once, with the understanding that his entire trip is one 
continuous circuit. There is no requirement that the distances be- 
tween cities be symmetric, that is, the distance from city i to city j 
need not be the same as that from j to i . The salesman 1 s objective is 

to select a sequence of cities in such a way as to minimize the total 

distance traveled. 

There are several algorithms which solve the traveling-salesman 

problem, but none offers the last word in terms of speed and ease of 

computation. The purpose of this investigation was to examine a newly 

1 

proposed solution procedure developed by Professor Harold Greenberg. 
His algorithm looked very promising in terms of hand computation on 
small problems, but this study focused on its desirability for 
handling larger problems, such as the forty-city case. The greatest 
effort went into programming the algorithm for the computer, then 
letting the computer serve as the primary tool of analysis. 



■^Greenberg, Harold, Professor, Department of Operations Analysis, 
Naval Postgraduate School, Monterey, California 
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A. DEFINITIONS 



The distance (cost) between city i and city j is denoted c— for 
i = l,...,n and j = l,...n. In general, c^j need not equal c__. The 
decision variable equals one if the salesman is to travel directly 

from city i to city j and equals zero, otherwise. Since traveling 
directly from city i to city i is not meaningful, x. . must always 
equal zero for all i. Most solution procedures meet this requirement 
by assigning arbitrarily large values to c^ for all i. 

Associated with each city is the number of steps required by a 
solution to arrive at that city. For example, if the salesman is to 
go from city A to B to C and back to A, he arrives at B in one step, C 
in two steps, and back to A in three steps. Unless otherwise stated, 
the following discussion takes city one to be the city of origin and 
city n to be the last city visited prior to returning to city one, 

A solution must provide a tour in order to be a traveling- salesman 
solution. A tour is defined as a sequence of transitions from city 
one through all other cities and back to city one, such that each city 
is entered exactly once and left exactly once, and it must take exactly 
n steps to return to city one. The latter clause guarantees there are 
no discontinuities in a tour. Deleting this clause yields the defini- 
tion of a subtour . It is possible for a candidate solution to contain 
as many as n/2 subtours and still meet the requirement that each city 
is entered and left exactly once. To illustrate the concept of tours 
and subtours, consider these solutions to the four city problem: city 

1 to 3 to 4 to 2 to 1 provides a tour; whereas city 1 to 2 to 1 and 
city 3 to 4 to 3 provides two subtours. 
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B . ALTERNATE FORMULATIONS 



The traveling-salesman problem is essentially a combinatorial 
problem and there is no one way to solve it. The most direct approach 
is to select a starting point, evaluate the distance (cost) over each 
of the (n-1)! possible tours, and then choose the tour(s) with minimum 
total cost. A dynamic-programming formulation reduces the number of 
tours one must consider but the tables required can quickly exceed 
computer storage capacity, Branch-and-bound techniques can also be 
employed but suffer from similar limitations as dynamic programming. 

A fourth possibility might be to construct a solution procedure using 
an algorithm for finding minimal-cost flows in a network. 

The traveling- salesman problem can be stated as an integer program- 
ming problem. Hadley j_ Ref, §J anc i Dantzig /"Ref. 1/ present essential- 
ly the same formulation, Dantzig lets x jk t = 1 or 0 according to 



whether or not the salesman travels from city i to city j on the t 

step, where i,j,t = l,...,n. Then, defining x. . , . = x. . , he 

ijl, 

states the problems as follows: 



th 



minimize 



^ — ' c x = z 
i>j,t ij ijt 



subject to 



T X --*. = I, x - , (j’ t= 

1 xjt k j , k , t+1 



and 



>_ x =1 (i = 1, . . . ,n) 

i , j , t 



and x ( , = 0 or 1 for all i,j,t. 

ijt ,J ’ 

The first set of constraints requires that if the salesman arrives at 

city j in t steps, he must leave that city on step t+1. The second 

sets requires that the salesman leave each city i once and only once. 

3 2 

Solving this system of n +1 variables and n + n equations with an 
integer programming algorithm yields a minimum-distance tour. 
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The integer programming approach has several shortcomings. The 
third subscript t on x greatly increases the number of variables, and 
the number of constraints in the formulation is very large. Further- 
more, integer programming algorithms normally require many more itera- 
tions than ordinary linear programming techniques » 

There is a general class of indirect solution procedures that 
takes advantage of a special feature of the traveling-salesman problem. 
If one temporarily ignores the fact that a solution must provide a 
tour, the salesman's task can be formulated as a simple assignment 
problem, perhaps the most readily soluble of all linear programs. 

The assignment formulation is achieved by letting x = 1 or 0 

iJ 

according to whether or not the salesman is to travel from city i to 
j. Let c be arbitarily large for all i. Then, 



minimize 




c x 

ij ij 



subject to 




(j + 1 






and x =1 (i = l,...,n) 

j ij 
2 

Notice that there are only n variables and 2n constraints. One 
constraint may be dropped because exactly 2n-l are linearly independent. 
Furthermore, since the assignment problem is a form of the transporta- 
tion problem, and since all restraining expressions equal one, any 
solution should yield 0 or 1 values for the x^ ^ 1 s /_Ref. 2 /, Thus, 
there is no need to add the integer constraint to the formulation 
above. Possible solution procedures include the simplex method, the 
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familiar transportation tableau, or an algorithm for finding minimal- 
cost flow in a network /_Ref. 4/. 

An assignment outcome does not necessarily provide a traveling- 
salesman solution, but it does satisfy the requirements that each 
city is entered exactly once and departed exactly once and that x^= 0 
for all i. Clearly, since the salesman 1 s problem includes the added 
constraint that any set of values for the x ij’ s must provide a tour, 
his minimum- distance solution must be at least as great as the minimum 
assignment distance. 

The so-called indirect procedures, alluded to above, typically use 
an assignment solution as a starting point. The optimal value of the 
assignment objective function provides a lower bound to any traveling- 
salesman solution. A non-optimal assignment solution may still provide 
a lower bound while being closer in magnitude to the salesman’s minimum 
total distance. In the sense that it is closer, the latter lower 
bound is said to be ’’better” than the former. In practice, however, 
it is difficult to find a value larger than the optimal assignment 
solution and prove that it is, in fact, a lower bound to the traveling 
salesman problem. 

Professor Greenberg's algorithm uses the optimal assignment solu- 
tion as its point of departure because it provides a guaranteed lower 
bound and because special use is made of the coefficient vectors in 
the final simplex tableau. After the solution of the initial assign- 
ment problem, the algorithm simply manipulates the coefficient 
vectors until a minimum-cost tour is achieved. This procedure is 
explained in detail in the next chapter. 
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W. L. Eastman also begins with the optimal assignment solution and 
if it should provide a tour, then the problem is solved. If, in the 
more likely event, it does not, he selects the subtour with the least 
number of variables and sets each equal to zero, one at a time. For 
each case, he resolves the original assignment problem but with the 
added constraint. In practice this can be achieved by setting 
c^j = CO corresponding to the required to be zero. If a tour is 
not obtained, he branches from that assignment problem with the 
smallest-valued objective function. Again the least-membered subtour 
is chosen to create new assignment problems, and so forth, until a 
tour is achieved. /Ref. 8/ 

J. D. C. Little and others have devised a branch-and-bound 
algorithm that starts with a lower bound, not to the traveling salesman 
problem, but to the assignment problem. Instead of obtaining the 
optimal solution to the assignment problem at each step, the smallest 
element is subtracted from each row of the cost matrix, and then the 
smallest element from each column of the result. The sum of the re- 
maining elements is said to be a valid lower bound on the optimal 
assignment solution. At each branching step two new problems are 

created corresponding to x . . = 0 and x. . = 1, The x. . chosen for each 

ij ij 

branching step is the one that yields as large a bound as possible 
when set equal to zero. Branching proceeds until a tour is obtained. 
Little ! s method normally requires many more branches than does Eastman* s, 
but there is less computation at each step. Some performance results 
are available on Little *s algorithm. The mean execution time, on an 
IBM 7090 computer for 100 thirty-city problems was 58,5 seconds and the 
average for five forty-city problems was 8,37 minutes. /Ref. 9/ 
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C. PLAN OF THE REPORT 



Chapter II presents a detailed description of the solution procedure 
suggested by Professor Greenberg. His proposed algorithm is complete 
in the sense that the criteria for proceeding from one step to the 
next are completely specified and guarantee an optimal solution. In 
some instances, however, it is not readily apparent how one can most 
efficiently meet these criteria. Chapter III deals with these ques- 
tions of technique as it describes how the algorithm was programmed 
for the computer. When there were alternate computational approaches 
possible, an attempt was made to find a good one. There is no claim 
that the particular techniques devised to meet the algorithm criteria 
are the ones Professor Greenberg, himself, would have chosen. 

Chapter IV outlines the salient results of computational experience 
with the computer program. The algorithm and the program are evaluated 
in the light of this experience. 
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I I . THE ALGORITHM 



An algorithm is defined by Webster as n a rule or procedure for 
solving a mathematical problem that frequently involves repetition of 
an operation. M This chapter seeks to present Professor Greenberg's 
algorithm for solving the traveling-salesman problem in the logical 
sequence of its development, concluding with a concise statement of 
the algorithm, itself, 

A. THE TRAVELING- SALESMAN PROBLEM 

A traveling salesman is faced with the problem of finding the 
shortest route (or the least-cost route) through n cities. Further- 
more, he must meet the following conditions: 1) he must enter each 

city once and only once, 2) he must depart from each city once and 
only once, and 3) his route must be a tour, that is, there must be no 
discontinuities in his route. There are (n-1)! possible solutions to 
this problem, since the salesman can go to any of n-1 cities on his 
first step, then to any of the remaining n-2 cities on his second 

th 

step, and so forth, until he returns to his starting point on the n 
step. Optimal solutions, of course, are found within this set of 
(n-1)! tours. 

Clearly, the above is merely an allegorical statement of a problem 
that has relevance to a variety of situations. For example, suppose 
fixed electronic components must be connected in a continuous circuit 
such that the least amount of wire is required. Disregard of condi- 
tions (1) and (2) might leave some components out of the circuit 
entirely or, on the other hand, cause short circuits. Violation of 
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condition (3) would result in two or more subcircuits, with current 
flowing in only one. 

B. THE ASSIGNMENT FORMULATION 

Modifying the traveling- salesman problem by disregarding the third 
requirement above, yields a problem which is readily adaptable to the 
assignment formulation, perhaps the most easily soluble of all program- 
ming problems. The assignment problem is usually stated in terms of 
assigning n men to n jobs in such a way that overall productivity is 
maximized, assuming that a productivity index can be associated with 
each man- job combination. There are n! possible solutions to this 
problem, since the first man can be assigned to any of the n jobs, the 
second man to any of the remaining n-1 jobs, and so forth. 

Dropping the requirement, then, that a solution must provide a 
tour, the modified salesman problem is to assign n cities to n cities 
(the same n cities) such that the sum of city-to-city distances is 
minimized. The possibility of assigning city i to itself is precluded 
by giving each c^ an arbitrarily large value. The possibility of 
other types of subtours appearing in a solution, however, cannot be 
restricted without destroying the assignment formulation. 

Since no city is ever assigned to itself in the modified salesman 
problem, the size of the solution set is smaller than n!, the usual 
size for assignment problems. The table below illustrates solution- 
set sizes for several selected problems: 
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Number 


Number of 


Number of 


Number of 


>f cities 


Assignment 


Modified Traveling- 


Solutions that 


(n) 


Solutions 


Salesman Solutions 


are Tours 




(n!) 


n! > (?)£• (n-1) ! 


(n-1)! 


2 


2 


1 


i 


3 


6 


2 


2 


4 


24 


9 


6 


5 


120 


44 


24 



A generalized expression for the third column above was not derived, 
but it is clear that the fraction of the time a modified salesman 
solution is also a tour is greater than (n-l)!/n! = 1/n. 

An optimal solution to the modified salesman problem provides a 
lower bound on the minimum-distance solution to the original problem. 
This is true, because the modified problem differs from the original 
only in that it has fewer constraints. If the optimal assignment 
solution to the modified problem provides a tour, then the original 
problem is solved; if not, further steps must be taken. 

To express the assignment formulation in the notation of Chapter I, 
let x^j = 1 or 0 according to whether or not the salesman is to travel 
from city i to city j; i,j = l,...,n # Let c. # be the distance between 
city i and j for all i and j. Set c ^ arbitrarily large for all i. 
Then, 



minimize 



— c . x 

j ij 



subject to 



^ x. . = 1 ( j = l, . . . ,n) 

i ij 



and 



x =1 (i=l , . . . ,n) 

ij 



and 



x. > 0 for all i , j 

x j 
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Alternately, this formulation can be expressed in matrix notation 
as follows: 

T 

minimize c x = z 

subject to Ax = b 

and x ^ 0 



where the definitions of c ,x,b, and A can be indicated most con- 
veniently by example. Consider the three-city problem: 

T 



x 



and A = 



1 

0 

0 

1 

0 

0 



c 12> « 


= 13’ C 21’ C 22’ C 23’ c 31» C 32’ 


C 33 ) 


X 11 




1 


x 12 




1 


x 13 




1 


X 21 




1 


X 22 


, b 


1 


x 23 




1 


X 31 




1 


x 32 




1 


x 33 




1 

, 



1 

0 

0 

0 

1 

0 



1 

0 

0 

0 

0 

1 



0 

1 

0 

1 

0 

0 



0 

1 

0 

0 

1 

0 



0 

1 

0 

0 

0 

1 



0 

0 

1 

1 

0 

0 



0 

0 

1 

0 

1 

0 



0 

0 

1 

0 

0 

1 



In general, c is 1 x n , x is n x 1, b is n x 1, and A is 2n x n 



th 



Notice that since the k un column of A corresponds to the °f x * 
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the i element and the (n+j) element of the k column of A equal 
one and all other elements in that column equal zero. 

The assignment problem has several important features that are 
listed below without proof . (see references 3 and 5 for proofs.) 

1. The assignment problem has a feasible solution, 

2. The rank of the matrix A is 2n-l, thus any row may be deleted 
as redundant, 

3. Any set of 2n-l linearly independent column vectors from A 
forms a basis B. 

4. The inverse of any basis B is composed entirely of zero, one 
and minus-one elements. 

5. Any basic feasible solution to Ax = b is degenerate. In 
fact, there are exactly n basic variables equal to one and exactly 
n-1 equal to zero, assuming all nonbasic variables are arbitrarily 
set equal to zero. 

C. THE SIMPLEX SOLUTION OF THE ASSIGNMENT FORMULATION 

There are alternate methods for solving the assignment problem but 
the following discussion assumes some variant of the simplex method is 
chosen. The final simplex tableau contains not only the optimal solu- 
tion to the assignment problem, but other interesting information, as 
well. Before proceeding, a few more notational conventions and 
definitions are stated below: 

1. Redefine the coefficient matrix A such that one row (any row) 
is deleted. 

2. Temporarily redefine the vector x such that the subscripts on 

2 

its elements run from 1 to n rather than 1,1 to n,n. 
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3. Make the same change in subscript notation for the elements 

* T 
of c . 

, th 2 

4. Denote the j column of A as a^ for all j = 1,..., n . 

5. Let B denote a basis in A, that is, any set of 2n-l column 
vectors/in A. 

6 . Let denote the vector of x j ' s corresponding to the a.'s in 
B. Call the elements of the values of the basic variables. Note, 



-1 



since x^ = B b, x^ is simply the vector obtained by summing each row 
of B, Requiring x^> 0, then, insures its elements are ones and zeros, 
exclusively. 

T 

7. Let c^ denote the vector of c j ' s corresponding to the elements 



in Xg. 



-1 



8 . Define y^ = B a . Note that the y . 1 s are composed entirely 
of zero, one, and minus-one elements. 

T 

9. Define the relative cost for x. as c. - c -c^y . 

n - * D - 



J 



J 



The simplex procedure involves finding a basis B such that x = 

B 



-1 



B b minimizes z. A feasible solution is recognized when the elements 
of xg are zeros and ones, exclusively, and an optimal solution, when 

w 2 

cj - 0 for all j = 1 , . . . ,n . 



Not 



ice that c^ = 0 always for all x j 1 s in x g because for 



such 



x j 1 s 3 c. = c b Yj • If Cj = 0 at optimality for a nonbasic x^ , then an 
alternate optimum can be obtained by bringing that Xj into the basis. 
In this case, the simplex methodology specifies which basic variable 
is to be replaced in order to maintain feasibility. 

To appreciate the impact of the relative costs on the value of z, 
consider the following relationship: 



^ = z + 21 c .x . 

j€R J J 
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where z is the new value of z and R is the index set of all nonbasic 
variables. This equation shows that if nonbasic variables are set 
equal to one, z is increased by the sum of their relative costs. 

Of course, if one or more nonbasic variables are set equal to one, 
some change must take place in the vector. This is true since an 
assignment solution requires exactly n variables equal to one, includ- 
ing the nonbasic variables. In order to measure such effects on x^., 
consider the following equality: 

x n = x - y x 

B B j£R 'j j 

A 

where x^ is the new set of values for x^. Normally x. = 0 for all j in 

R, but if one or more such nonbasic variables are set equal to one, the 

relation above changes the values of the basic variables. Setting one 

or more x. f s not in x_ equal to one is said to be constraint feasible 

as long as the elements of x are zeros and ones, exclusively. 

B 



D. STATEMENT OF THE ALGORITHM 



The algorithm under study begins with the optimal assignment 
solution to the modified traveling- salesman problem. If the optimal 
assignment solution provides a tour, then the original problem is 
solved. Otherwise, the algorithm begins to introduce combinations of 
nonbasic variables at the one level until a tour is achieved. 
Complexities arise in the measures taken to assure that 
combinations are considered in ascending order of their relative costs. 
If these steps are not taken, there is a possibility that nonbasic 
variables not yet tried yield lower-cost tours. 

In the light of the background developed in the previous sections, 
the algorithm is as follows: 
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1. Drop the constraint that a traveling-salesman solution must 
provide a tour, 

2. Solve this modified problem as an assignment problem using a 
variant of the simplex method. 

3. Check the optimal assignment solution for a tour. If a tour 
is provided, stop. Otherwise, continue. 

4. Drop simplex tableau vectors corresponding to the original x^ 

variables. Drop tableau vectors corresponding to basic variables. 

2 2 
There are n -3n+l vectors remaining. Let m = n -3n+l. 

5. Rank tableau vectors in ascending order of their respective 
relative costs. 

6. Let x. = x.,y. = y., and c. = c. where i = 1 if c. is first 

J i ,J j j i J 

in the ranking, i = 2 if Cj is second, and so forth, through i = m if 
Cj is last in the ranking. (Retain information required to reconvert 
subscripts as needed when checking for tours.) 

7. Let i = 0 # 

8. Increment i by one. 

9. If x. = 1 is not constraint feasible and i = 1, return to 
step 8. If x- = 1 is not constraint feasible and i>l, go to step 11. 

If x- = 1 is constraint feasible, continue. 

10. If x. = 1 and x_, = x -y provide a tour, go to step 21. Other- 

1 D D i 

wise, return to step 8 if i = 1; continue if i >1. 

11. Generate all combinations of x, = 1 with x = 1, x =1,..., 

i 1*2 

x^_^ = 1. (Note this is construed to mean x^ = 1 appears in each 
combination along with one or more of the other variables indicated set 
equal to one. There are 2^ ^ - 1 such combinations.) 
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12. Disregard all combinations with relative costs greater than 

c i+1 . (Define - ^ ,) 

13. Disregard all remaining combinations that are not constraint 
feasible. 

14. Check remaining combinations, in ascending order of their 
relative costs, for tours. As soon as a combination is found to pro- 
vide a tour, go to step 21. If no combination yields a tour, continue 

15. If i = 2 and/or c^ + ^ = c^, return to step 8. Otherwise^ 
continue. 

16. Let k = i. 

17. Decrement k by one. 

18. Generate all combinations of = 1 with = 1, ^ = 1,..., 
\-i = 1 - 

19. Same as steps 13, 14, and 15. 

20. If k > 2, return to step 17. If k = 2, return to step 8. 

21. If coming from step 19, let i = k. Otherwise / do not change i 

22. Compute total distance for the tour found. Letting K be the 

index set of nonbasic variables included in the solution, the total 

distance for this tour is z = z c, . 

k4K k 

23. If i ^ 2, go to step 32. Otherwise, continue. 

24. If all combinations of nonbasic variables preceding x^ have 
already been checked against the current value of z-z as an upper 
bound, go to step 32. Otherwise, continue. 

25. Decrement i by one. 

26. Generate all combinations of x. = 1 with x, = 1, = 1, . . . , 

i 1 2 

x . = 1 . 

i-l 

27. Disregard all combinations with relative costs greater than 

A 

z-z. 
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28. Disregard all remaining combinations that are not constraint 
feasible. 

29. Check remaining combinations, in ascending order of their 
relative costs, for tours. As soon as a combination is found to pro- 
vide a tour,- go to step 30. If no combination yields a tour, go to 
step 31. 

30. Update the value of z for the latest tour found and continue, 

31. If i y 2, return to step 25. Otherwise, continue. 

32. The latest tour found and its associated z is an optimal 
traveling-salesman solution. Stop. 

Since eventually all combinations of nonbasic variables are con- 
sidered, the algorithm must eventually find a tour. Suppose a tour 
is first found using x. = 1 in some combination with the preceding 
nonbasic variables* This tour’s relative cost is less than or equal 
to c^_^ because only combinations meeting this criterion were considered. 
Thus, there is no need to try combinations beyond x^, There is a need, 
however, to reconsider combinations prior to x^. For example, x^ = 1 
in combination with x^ = 1 may have been rejected earlier because its 
relative cost is greater than c^, but it may be less than z-z. Thus, 
proceeding backward through the nonbasic variables using the ciirrent 
value of z-z as an upper bound guarantees the minimum 4 will be found. 
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III. THE COMPUTER PROGRAM 



The algorithm was programmed for the IBM System/360 Model 67 
Operation System using FORTRAN IV (H level) programming language 
/Ref. 7/ . A listing of the complete program is provided in Appendix 
A, 



Several routines used in the program were not developed as part 

T 

of this study but were modified for use here as needed. Subroutines 
ASSIGN, LIP, ELEM, and PIVOT were supplied by Professor Greenberg and 
are used to solve the initial assignment problem. Subroutine RAND is 
a simple routine for generating city-to-city distances, uniformly 
distributed from 0 to 100. The distance from each city to itself is 
arbitrarily set equal to 999. Subroutine TIMEIT is used to measure 
algorithm execution times in seconds. Part of TIMEIT appears in 
assembly language. 

The following list gives a brief description of the various 
routines developed especially for this study. 

1. The Main Program follows the sequence of algorithm steps as 
stated in Chapter II. Provision is made to either read in a cost 
matrix or to generate one randomly. Other provisions are made to 
abort if something unexpected goes wrong and to begin execution on 
the next problem. 

2. Subroutine PRINTC writes out the cost matrix. Also, alternate 
storage space is provided for the costs because they get overlayed 

in the assignment routines by the relative costs. 

3. Subroutine PRINTA writes the optimal assignment solution in 
tableau form. If the problem is larger than seven cities they y^ 
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column vectors are not printed because there are too many to fit on a 
page. Subroutine PRINTA is also used again to print the optimal assign- 
ment solution after artificial variables remaining in the basis have 
been replaced, 

4. Subroutine SOLCHK checks constraint feasible solutions for 
tours. If a tour is found, it t is printed. 

5. Subroutine RANKC ranks the relative costs in ascending order of 
magnitude. Actually, the costs are not moved around in storage, but 
second-order subscripts are assigned to indicate the desired ranking. 

6, Subroutine REPART was devised to replace any artificial 
variables left in the assignment solution and to select a redundant 
row for deletion from the final tableau. 

7, Subroutine LINCOM is the most important subroutine in terms 
of the algorithm at hand. This subroutine generates the combinations 
of the nonbasic variables, compares the relative cost of each combi- 
nation with an upper bound, and checks for constraint feasibility. 

A. DEFINITIONS 

Three common blocks (storage areas) are used in the program. 

Common Block ASGBLK is used exclusively with the four assignment 
routines. Common Block SHARE is used to pass information from the 
assignment routines and from Subroutine RAND to the traveling- sales- 
man algorithm routines. Common Block TSBLK is used to pass values 
among the algorithm routines. For easy reference, this section 
defines the variables in common in the order of their appearance. 
Variables in ASGBLK are not listed because they are only needed 
within the four assignment routines which are not part of this 
study. 
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1 . Common Block TSBLK 



a, COST is the sum of the relative costs included in a 
particular combination of nonbasic variables, each taken at the one 
level. 

b. CSTLC(1600) is initially used as alternate storage for 

the city-to-city costs. After artificial variables have been replaced, 
it is used to store relative costs of combinations in ascending order 
of magnitude. 

c„ NLC(IOOO) is used to store a combination identifying 
number corresponding to each value in CSTLC(1600). With these numbers, 
combinations can be regenerated as needed. 

d. IC(1600) is used to store the subscripts of the variables 
associated with the ranked, relative costs. For example, the first 
value in IC(1600) is the subscript of the nonbasic variable with the 
smallest relative cost. 

e. KC(40) is used in Subroutine PRINTC to write out the city- 
to-city costs as integers. 

f. KB(120) is used to store the new values of the basic 
variables when one or more nonbasic variables are set equal to one. 

The one-values for the nonbasic variables are augmented to this vector 
starting after the last basic-variable valua (Also, see description 
of IB(120) below.) 

g. IAX(40) is used in Subroutine PRINTA to write out a row 
of the final assignment tableau. 

h. IROW(40) and ICOL(40) are used in Subroutine SOLCHK to 
store the first and second subscripts, respectively, of variables after 
these have been reconverted from single-subscript form. 



24 



i. LC(40) is used to store a binary number up to 40 digits in 

length. If the digit is a: one, then the ranked, nonbasic 

variable is to be included in the particular combination at the one 

level. If the k digit is zero-, then the k t ^ 1 ranked, nonbasic variable 

40 

is not to be included, LC(40) can be used to represent up to 2 - 1 

different combinations, one at a time. 

j. KZ is the truncated value of Z in Subroutine PRINTA. In 
Subroutine SOLCHK, KZ is the truncated value of XZ = Z + COST. If the 
city-to-city costs are integers, KZ = XZ always. 

k. NC = N-3MM+1 is the number of nonbasic variables, excluding 
x ii var i a kl es an ^ artificials, 

l. J is used in the Main Program to index the progression 
through the ranked, nonbasic variables. 

m. ICJ = IC(J). 

n. Ll - LEG(1, ICJ ) and L2 = LEG(2,ICJ). 

o. Jl = J-l. 

p. CST is the current upper bound against which combinations 
are compared. Any combination of nonbasic variables at the one level 
with a relative cost greater than CST is disregarded. 

q. XZ = Z+COST. Subroutine SOLCHK truncates this value when 
writing it out. (See definition (j) above.) 

2 , Common Block SHARE 

a. C(1600) is used initially to store the city-to-city costs. 
After the assignment solution, the relative costs are stored here. 

(Also, see description of CSTLC (1600) above.) 

b. B(80,80) is the inverse of the optimal assignment basis. 
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c. H(80) is used to store the values of the basic variables 
in the optimal assignment solution. (Also, see description of KB(120) 
above. ) 



d. Z is the optimal assignment total cost. 

e. LEG(2,1600) is a matrix of values used to convert single- 
subscripted variables back to their original double-subscripted form. 
For example, x, = x. . if i = LEG(l,k) and j = LEG(2,k) - MM where MM 

x ij 

is the number of cities. Furthermore, LEG provides a convenient device 

- 1 f v, 

for obtaining the product B a^ = y^.. It is known that the i and 

th 

(j+MM) elements of a^ are ones, if i and j are defined as above, and 
all other elements of a^ are zeros. Thus, rather than multiplying 
B 1 times a , simply add the i^* 1 and (j+MM) ^ elements of each row of 

K 



B to produce the vector y^. 

f. IB(120) is used to store the subscripts of the optimal 
assignment basic variables. Negative subscripts indicate artificial 
variables. Starting just after the last basic-variable subscript, 
nonbasic-variable subscripts are appended to this vector for each 
nonbasic variable set equal to one in a particular combination. (Also, 
see description of KB(120) above.) 

g. MM is the size of the problem, that is, the number of 



cities. 



h. M is the number of constraining equations in the assignment 
formulation. Initially, M = 2 (MM) , but in Subroutine REPART it is 
reduced by one. 

2 

i. N = (MM) is the total number of variables, excluding 
artificials . 
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B. DESCRIPTION OF SELECTED SUBROUTINES 

The Main Program was written to conform closely to the algorithm 
as presented in Chapter II and can be followed readily from the listing 
in Appendix A. Further amplification of Subroutines LINCOM, SOLCHK, 
and REPART is considered necessary, however, in order to clarify the 
workings of the overall program. Subroutine REPART was devised only 
because- the assignment routines used leave artificials in the basis 
and do not drop a redundant equation. Subroutines SOLCHK and LINCOM 
represent techniques created to fulfill algorithm requirements and, as 
such, are the most important products of this study. 

'* 1 • Subroutine REPART 

Subroutine REPART scans the subscripts of the basic variables 
returned from the assignment routines. As soon as a negative sub- 
script is located, its corresponding artificial variable is replaced 
in the basis by^.the first, ranked, nonbasic variable that has a nonzero 
element in the corresponding row of its vector. Using this nonzero 
element as a pivot element, standard row transformations are performed. 
The result is that the nonbasic variable is brought into the basis and 
the artificial is dropped. The optimal value of the objective function 
does not change because the artificial variable was necessarily at the 
zero level implying its replacement must also be at the zero level. 

After replacing an artificial, relative costs for the remaining non- 
basic variables are recomputed and these variables are then ranked in 
ascending order of their new relative costs. The above process is 
repeated until all artificials have been removed from the basis. 

It is known that one of the assignment- tableau rows is redundant. 
If, while trying to replace an artificial, a nonzero pivot element 
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cannot be found, that row is earmarked for deletion. Otherwise, after 
all artificials have been replaced, each row is scanned for all zero 
elements. If such a row is found, it is dropped,, Otherwise, the last 
tableau row is arbitrarily selected for deletion. If other than the 
last row is dropped, remaining rows are shifted upward in the tableau 
to fill the gap. 

Normally, control is returned to the statement following the 
calling statement in the Main Program. If for any reason, however, 
more than one row has gone to all zeros, control is returned to state- 
ment 10 in Main and execution on the next problem is started, 

2, Subroutine SOLCHK 

Subroutine SOLCHK examines an assignment solution to see if it 
provides a tour. Its approach is to count the number of steps required 
to return to the city of origin. If this number is less than the total 
number of cities MM, then subtours exists and control is returned to 
the statement following the calling statement in Main* If the number 
of steps required is exactly equal to MM, the solution is a tour. In 
the latter event, the solution is printed with its associated total 
distance and control is returned to the statement in Main indicated in 
the parenthesis of the calling statement* 

To illustrate the bookkeeping involved, consider the two fol- 
lowing four-city assignment solutions: 



X 

00 

II 


x 2 ,4 = 


x 3 “ 


X 1 , 3 = 


x 14 : 


= X 4,2 


II 

ON 

X 


x 3, 1 = 



Steps 

1 



2 
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Steps 

1 



3 



x 



2,4 



1 



2 



x, 



'2 



x l,2 = 1 



4 



Both are assignment solutions, but only the second provides a tour. 

Notice that the routine will keep scanning the KB vector until 

th 

the (MM) one-value is found and then scanning is terminated. This 
approach is justified because only constraint feasible solutions are 
sent to SOLCHK and these solutions necessarily have exactly MM 
variables equal to one, including any nonbasic variables set equal 
to one. 

3. Subroutine LINCOM 

The essence of the algorithm lies in altering the optimal 
assignment solution by setting selected nonbasic variables equal to 
one until a tour is achieved. Subroutine LINCOM was devised to generate 
all the required combinations of nonbasic variables. An attempt was 
made to program LINCOM so that the number of combinations necessary to 
consider be held to a minimum. The objective was to generate them in 
ascending order of total relative costs so that as soon as a combina- 
tion^ cost exceeded the current upper bound, no further combinations 
need be considered. This objective was only partially achieved. 

The most difficult problem lay in the fact that there is no 
way to completely predetermine the ranking of all possible combinations. 
Consider the case where c-^ ^ ^2 ~ • •• C 5~ ^6 anc * x 5 = ^ ^ ounc * not 
constraint feasible. Before testing Xg = 1, all combinations of 
x^ = 1 with the preceding nonbasic variables must be tried because some 
may have total relative costs less than or equal to c^. 
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One way to test the combinations would be to try x^ = 1, Xg = 1, 
and =. 1 singly with = 1, then in pairs, then triplets, and 
finally all four, A better approach, however, is to set up the four- 
digit binary number 0000 and then simply increment it by one, fifteen 
times. Number the digit positions from one through four, from right to 
left. Then, if a one appears in the position where j = 1,,..,4, 
let x_. = 1 appear in the combination. The possible combinations can 
now be represented as shown below: 

Identifying 



Number 


Combination 


i 


0001 


2 


0010 


3 


0011 


4 


0100 


5 


0101 


6 


0110 


7 


0111 


8 


1000 


9 


1001 


10 


1010 


11 


1011 


12 


1100 


13 


1101 


14 


1110 


15 


1111 



Combination number 9, for example, represents x,. = 1 with x^ = 1 and 
x-l = 1 for a total relative cost of c 5 + 5^ + c^. 

Notice that the following generalizations can be made about 
the respective total relative costs of the combinations above: 





Ml 

I — 1 


#2 




#3, 




#4 — 


#5 




#6 — #7, 




Ml 

00 


#9 




#10 ^ #11, 


and 


Ml 

CNl 
!— < 


#13 


— #14 ^ #15, 
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Furthermore, #3 — #5 ^ #9, but #3 1 s relationships with #4 and #8 

cannot be stated, in general. Also, #7 ^ #11 ~ #13, but #7 1 s 

relationships with #8, #9, #10, and #12 are unknown. Finally, #11^: 
#13, but there is no way to tell in advance which is smaller, #11 or 
# 12 . 



While there are always many organizational difficulties 

associated with writing and debugging a large program, the situation 

discussed above presented one of the most challenging theoretical 

problems to resolve. The scheme finally chosen is stated below: 

If a combination's relative cost is greater than the current upper 
bound and the combination's identifying number is equal to 

i 

1) 2 , where i is any nonnegative integer, terminate genera- 
tion of further combinations (because all their relative 
costs will exceed the current upper bound), or 

2) 2 1 + 2^, where i and j are any nonnegative integers such 
that j < i, skip to combination number 2 1_ *~1 and continue 
generating, or 

3) otherwise, continue generating. 

In most cases, the above procedure greatly reduces the number of 

combinations necessary to generate. It could be expanded by saying 

i j k 

that if the combination number is equal to 2+2+2 , where i, j, and 
k are nonnegative integers such that k<j<i, skip to number 2^+2 
and so forth. Except for very large problems, however, such added 
checks would probably consume more computer time than they would save. 

Another facet of the problem of trying to generalize the 
relationships among relative costs is presented when a combination is 
found constraint feasible. To illustrate the situation, suppose 
combination number 11 is found to provide a relative cost less than the 
current upper bound and is also found constraint feasible. If #11 is 
sent to SOLCHK at this point and found to provide a tour, there can be 
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no certainty that it is a least-cost tour among the current set of all 
possible combinations. This is because #12, as yet unchecked, may 
provide an even smaller relative cost, may be constraint feasible, and 
may also provide a tour. 

The technique devised to resolve the above difficulty was to 
store the identifying number and relative cost of each combination 
found constraint feasible. With each addition to storage, the numbers 
and relative costs are reranked in ascending order of relative costs. 
Finally, when all necessary combinations have been generated and 
checked, the constraint feasible solutions are sent to SOLCHK in 
ascending order of relative costs, thereby insuring the first tour 
found has a cost less than or equal to any tour the remaining solutions 
in storage might provide. As soon as a tour is found, control is 
returned to the first address indicated in the parenthesis of the 
calling statement in Main, If none of the solutions in storage provides 
a tour, control is returned to the statement in Main following the 
calling statement. 

Finally, a word about the regeneration of combinations required 
by the above procedure. The decimal identifying number in storage is 
converted to binary through the standard technique of successive 
divisions by 2. Remainders become the binary digits from right to left. 
The binary number, then, is the original representation of the combina- 
tion represented temporarily in storage by the decimal identifying 
number. Using the binary number, the combination of interest is re- 
constructed and sent to SOLCHK along with its associated relative 
cost. 
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C. LIMITATIONS OF THE PROGRAM 



The storage areas of the program are set up to deal with problems 

smaller than 41 cities, but these areas can be expanded easily if 

necessary. The four subroutines used to provide an optimal assignment 

solution are programmed to terminate execution in the unlikely event 

that more than 299 iterations are required. The only other significant 

limitations of the program are found in Subroutine LINCOM, but one of 

these presents a serious restriction. 

In Subroutine LINCOM, up to 1000 constraint feasible solutions can 

be stored at one time, but this does not appear to present an important 

limitation. In any case, storage space for CSTLC and NLC in Common 

Block TSBLK can be increased considerably if necessary. The most 

serious restriction and the one hardest to correct is the fact that 

Subroutine LINCOM can generate all the possible combinations for no 

more than 32 nonbasic variables. This means that if a tour is not 

r d 

found upon setting the 33 ranked, nonbasic variable equal to one, 

then execution on this problem must terminate. The reason for this 

situation is that each combination must have a unique identifying number 

associated with it in LINCOM and the largest integer number that can 

31 32 

be expressed on the IBM 360 is 2 -1. There are 2 -1 possible 

r d 

combinations of the 33 variable with the 32 preceding it, hence, 
the problem. As an example of the gravity of the restriction, note 
that a forty-city problem presents 1481 nonbasic variables for 
consideration. 

The perplexity above is not insurmountable, but it was not re- 
solved in this study. The possibility of representing the identifying 
number in floating-point mode is precluded because a form of the 
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number is used as a DO Loop parameter. (See the variable KT in LINCOM) 

It might be possible to represent the number as an element of a vector. 

31 

For example, each of the first 2 -1 numbers might be prefixed with a 

31 

1, the second set of 2 -1 numbers with a 2, and so forth. The DO 

Loop using the number as a parameter could, in like manner, be 
expanded into a double DO Loop, 
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IV. COMPUTATIONAL EXPERIENCE 



Unfortunately, time constraints on this study did not permit 
extensive testing of the computer program in its final form. Enough 
was accomplished, however, to indicate several important conclusions. 

Appendix B presents a partial facsimile of the program output for 
the same ten-city problem that is discussed in References 8 and 9. 
Eastman* s algorithm required the solution of eleven assignment prob- 
lems before the optimal tour was obtained. Little* s method required 
31 branches to arrive at the desired solution. Appendix B indicates 
the programmed algorithm under study took only 1.50 seconds after the 
initial assignment solution to obtain the optimal tour. The first 
tour was found after trying combinations of the first 23 ranked, non- 
basic variables and this tour turned out to be the optimal one. 
City-to-city costs were generated by Subroutine RAND for 17 

twenty-city problems. Fourteen of these problems were solved by the 
* 

*■ / 

program in a meantime of 25 seconds and after trying an average of 15 
nonbasic variables for each problem. Execution was terminated on the 
other three problems because a tour had not been found within the first 
33 nonbasic variables. 

Ten randomly generated thirty-city problems were attempted and 
nine were solved in an average time of 56 seconds. The mean number 
of nonbasic variables required for testing in these nine problems was 

16. 

Finally, Subroutine RAND was used to create three forty-city 
problems. The first ran .07 seconds after the assignment solution to 
discover this was also an optimal tour. The second problem ran 549 
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seconds before terminating after trial of the 33 nonbasic variable 

failed to provide a tour* The third problem was solved in 187 

seconds after trying only four nonbasic variables* 

In view of the computational experience cited above, one can 

conclude that the algorithm works and that the program works with one 

exception. Of the 31 problems tried, five could not be solved because 

r d 

tours had not been found upon setting the 33 nonbasic variables 

equal to one* If the program is to be of general use, this somewhat 

arbitrary restriction in Subroutine LINCOM must be removed. The 

limitation could be withdrawn if a way is found to obviate the 

identifying numbers presently associated with each possible combina- 
31 

tion. Since 2 -1 is the largest integer that can be represented on 

nd 

the IBM System/360, none beyond the 32 n ranked, nonbasic variable can 
be tried in combination with its preceding variables. 

Several steps can be taken to increase the execution speed of the 
program. All write statements, except those in Subroutine SOLCHK and 
the various termination messages throughout the program, can be replaced 
with CONTINUE statements. In subroutine REPART there is no particular 
reason why relative costs are recomputed and reranked after each 
artificial replacement. These two chores could be postponed until all 
artificials have been replaced. Since a replacement variable enters 
the basis at the zero level, it is not necessary that its relative 
cost be the least possible. 

In conclusion, a possible improvement to the algorithm, itself, is 
suggested. The algorithm is very efficient if a tour is found within 
the first few nonbasic variables tried, otherwise much time is consumed 
rechecking combinations already considered against a new upper bound, 



36 



and perhaps again against an even newer upper bound, and so forth. 
This phenomenon can be observed in the example output displayed in 
Appendix B. The algorithm might be improved as follows: 1) obtain 

a relatively low-cost tour by linking together, on some reasonable 
fashion, the subtours appearing in the optimal assignment solution, 

2) compute z for this tour, 3) using z-z as an upper bound, begin 
with the last, ranked, nonbasic variable whose relative cost is less 
than z-z and proceed backward just as in the original algorithm after 
the first tour is found. In this manner, it is not necessary to 
proceed forward through the ranked, nonbasic variables looking for 
that first tour, and it is necessary to work backwards through the 
variables only once. 
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APPENDIX A 



THE PROGRAM LISTING 



C PLACE GREEN JOB CARO HERE. 

// EXEC FORTCALG, REGION. G0=160K, TIME. G0=10 
//FORT. SYSPR INT DO SPACE=(CYL» (4,1)) 

//FORT . SYS IN DD * 

C A TRAVELING-SALESMAN ALGORITHM 

CCMMON /TSBLK/ CO ST, C STLC ( 1600 ) , NLC ( 1300 ) , 

1 IC( 1600) ,KC(40) ,KB< 120) , I AX (40) , 

2 IR0W(40 ) , IC0L(40 ) ,LC(40) , 

3 KZ ,NC, J, ICJfLl ,L2, J1 ,CST,XZ 
CCMMON /SHARE/ C(1600),B( 80» 80), H( 80), Z, 

1 LEG(2»1600) ,IB(120),MM,M,N 

C MM IS THE NUMBER OF CITIES. 

C NRAND=1 IF COSTS ARE TO BE RANDOMLY GENERATED. 

10 RF AD ( 5 »50C ,END = 90 ) MM , NRAND 
500 FORMAT (21 2) 

M=2*MM 

N=MM*MM 

IF(NRAND.EQ.l) CALL RANDU20) 

RE AD (5, 510) (C(J) ,J=1,N) 

510 FORMAT ( 1 6F 5. 0 ) 

20 CALL PRINT C 

C OBTAIN OPTIMAL ASSIGNMENT SOLUTION. 

CALL ASSIGNU10) 

CALL PRINTA 
C START CLOCK. 

CALL T IME I T( 0 »TI ME ) 

C SEE IF OPTIMAL ASSIGNMENT SOLUTION IS A TRAVELING- 
C SALESMAN SOLUTION. IF YES GO TO 80, IF NO CONTINUE. 
COST=C. 

CALL SOLCHK( C80 ) 

C RANK NONBASIC VECTORS IN FINAL TABLEAU IN ASCENDING 
C ORDER OF THEIR ASSOCIATED RELATIVE COSTS. 

CALL RANKC 

C REPLACE ARTIFICIALS REMAINING IN OPTIMAL ASSIGNMENT 
C SOLUTION AND DROP A REDUNDANT EQUATION FROM THE TABLEAU. 
C NOTE THAT M IS DECREMENTED BY 1 IN SUBROUTINE REPART. 
CALL REPART(CIO) 

CALL PRINTA 

C( IC ( NC+ 1 ) )=1QOOOOO. 

C STHLD=C . 

C NC=MM*MM— 3*MM+1 IS THE TOTAL NUMBER OF NONBASIC 
C VARIABLES OF INTEREST. 

C PROCEED FORWARD THROUGH THE RANKED, NONBASIC VARIABLES. 
DO 5C L=1,NC 
J = L 

I CJ= IC ( J ) 

KCST=C ( I C J ) 

WRITE(6,61C) IC J , KOST 

610 FORMAT (//1CX, ' TRY X ( * , I 2 , * ) = 1 REL COST = ',13) 

L 1 = LEG ( 1 , ICJ ) 

L2=L EG ( 2 , I C J ) 

I B ( M+ 1 )=ICJ 
KB ( M+ 1 ) = 1 
COST=C( ICJ ) 

DC 30 1 = 1, M 

XHOLD=H( I ) -B ( I , L 1 ) -B( I ,L2 ) 

K B ( I ) = XHOL D+S I GN ( .5,XH0LD) 

C IF X ( IC J ) =1 IS NOT CONSTRAINT FEASIBLE, GO TO 40. 

IF(KB( I ) .NE.O. AND.KB( D.NE.l) GO TO 40 
30 CONTINUE 

CALL S0LCHM&58) 

40 IF(J.EQ.l) GO TO 50 

C TRY ALL POSSIBLE COMBINATIONS OF X ( I C J ) = 1 WITH THE 
C PRECEDING NONBASIC VARIABLES. DISREGARD ANY COMBINATION 
C WITH RELATIVE COST GREATER THAN CST=C ( I C ( J+ 1 ) ) . 

J1=J-1 

WPITE(6,620) J1 

620 FORMAT ( 10X ,* TRY THE (2***, 12,') - 1 LN COM') 
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CST=C ( IC( J + L > ) 

C IF A COMBINATION PROVIDES A TOUR GO TO 58, 

C OTHERWISE CONTINUE. 

CALL L I NCOM ( £ 5 8 » £ 1 0 ) 

C IF C ST IS EQUAL TO THE PREVIOUS CST GO TO 50. 

IFtCST .EQ.C( I C J > ) GO TO 50 
IF(J.EQ.2) GO TO 50 

C PROCEEDING BACKWARD THROUGH THE RANKED, NONBASIC 
C VARIABLES, RECHECK OLD COMBINATIONS AGAINST THE NEW CST. 
C STHLD=CST 
K J 1= J— 1 
DO 45 K=2, KJ 1 
J = K J 1-K+ 2 
J1=J-1 
I C J= I C ( J ) 

L1=LEG(1,ICJ) 

L 2=L EG ( 2 » I CJ ) 

I B ( M+ 1 ) = IC J 

K B ( M+ 1 ) = 1 

WR IT E( 6, 630 ) ICJ 

630 FORMAT (//1CX, • RETURN TO X < • , I 2 , * I AND') 

WRITEt6,62C) J1 

C IF A COMBINATION PROVIDES A TOUR GO TO 60, 

C OTHERWISE- CONTINUE. 

CALL L INCOMt £60 » £ 1C ) 

45 CONTINUE 
50 CONTINUE 

C THE PROGRAM SHOULD NEVER REACH THIS POINT BECAUSE 
C THERE ARE (MM-1): POSSIBLE TOURS. 

CALL TIMEITI-l .TIME) 

WRITE! 6 , 64C ) 

640 FORMAT (//10X, • A TRAVELING SALESMAN SOLUTION*, 

1 • DOES NOT EXIST* ) 

WR IT E( 6, 65C ) TIME 

650 FORMAT (//IOX ,* ALGORITHM TIME =*,F10.4,' SECONDS') 

GO TO 1C 

C IF CLC COMBINATIONS HAVE ALREADY BEEN CHECKED AGAINST 
C AN OLD CST .LE. THE RELATIVE COST OF THE TOUR JUST 
C FOUND, GO TO 80. 

58 IFtCST. LE.CSTHLD) GO TO 80 
60 IE(J.LE.2) GO TO 80 

C PROCEEDING BACKWARD THROUGH THE RANKED, NONBASIC 
C VARIABLES, RECHECK OLD COMBINATIONS AGAINST THE NEW CST 
C WHICH IS THE RELATIVE COST OF THE LATEST TOUR FOUND. 

K J 1= J- 1 
DO 70 K=2 , K J 1 
J = KJ l-K+2 
J1=J-1 
I C J= IC ( J ) 

L 1 = L EG ( 1, ICJ ) 

L2=LEG(2, ICJ ) 

I B ( M+l ) = I C J 
K B ( M + l ) = 1 
WRITE(6,630) ICJ 
WR I TE ( 6 , 62 0 ) J1 

C WHETHER OR NOT A COMBINATION IS FOUND TO PROVIDE 
C A TOUR, CONTINUE. 

CALL LINCOMt £70,810) 

70 CONTINUE 
C STOP CLOCK. 

80 CALL TIME ITt -1 ,TI ME) 

WRITE(6,66C) KZ 

660 FORMAT t//10X,*** Z =',I5,' IS OPTIMAL **•) 

WRITEt 6, 650 TIME 
GO TO 1C 
DO STOP 
END 
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SUBROUTINE PRINTC 

COMMON /TSBLK/ CO ST, C S TLC ! 1 600 > , NLC ( 1000 ) , 

1 IC( 1600) ,KC<40) ,KBl 120 ) , I AX !40 ) , 

2 IROW(40) ,ICOL(40) ,LC(40) , 

3 KZ » NC » J » ICJ,L1,L2,J1,CST, XZ 
COMMON /SHARE/ C(1600),B( 80 , 80, H( 80), Z, 

1 LEG( 2, 1600) ,IB( 120) ,MM,M,N 

WRITE(6,6C0) 

6 ")0 FORMAT!’ 1' ,11X,'C0ST MATRIX') 

WRITE! 6,610) (I, 1=1, MM) 

610 FORMAT (//11X, 4013 ) 

WRITE! 6,620) 

620 FORMAT ( ’ ' ) 

1 = 0 

DO 20 K=1 , N» MM 
L=K+MM-1 
J = 0 

DO 10 LL=K,L 
J = J+1 

KC!J)=C!LL) 

10 CONTINUE 
1=1 + 1 

WRITE16.630) I , ( KC ( J ) , J=1 , MM ) 

630 FORMAT! 8X , I 2 , 1 X , 40 1 3 ) 

20 CONTINUE 

DO 30 J=1 , N 
30 CSTLC! J)=C(J ) 

RETURN 

END 



SUBROUTINE PRINTA 

COMMON /TSBLK/ COST , CSTLC ! 1600 ), NLC ! 1 000 ) , 

1 IC! 160 0) ,KC140) ,KB! 120) , I AX! 40) , 

2 IROW!40),ICOL!40),LC!40), 

3 KZ,NC,J,ICJ,L1,L2,J1,CST,XZ 
COMMON /SHARE/ C!1600),Bi 80, 80). Hi 80), Z, 

1 LEG<2, 1600 ) , IB! 120 ) ,MM,M,N 

WRITE! 6,600) 

600 FORMAT i • 1 ', 10X , 'OPTIMAL ASSIGNMENT SOLUTION') 

IF <M.EQ. !2*MM) ) GO TO 10 
WRITE! 6,610) 

610 FORMAT i 19X » ' (MODIFIED) ' ) 

GO TO 60 
10 K Z = Z 

DO 20 1=1, M 
20 KB ( I ) =H ! I » +. 5 
NC=0 

MM1=MM+1 
DO 50 J= 1 , N 
DO 3C 1=1, M 

IF(J.EO.IBd) ) GO TO 50 
30 CONTINUE 

DO 40 K=1,N,MM1 
IF(J.EQ.K) GO TO 50 
40 CONTINUE 
NC=NC+1 
I C ( NC ) = J 
50 CONTINUE 

60 I F ( MM. LE . 7 ) GO TO 80 
WRITE! 6,620) 

620 FORMAT (//1CX, • BASI S B'//) 

DO 70 1=1, M 

WRITE! 6,630 IB! I ) , KB! I ) 

630 FORMAT (1CX, 14, 2X, 12) 

70 CONTINUE 

WR ITE ( 6, 64C ) KZ 
640 F0RMAT(/15X* *Z = ', 14) 
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RETURN 

80 WRITE(6,650) 

650 FORMAT !//22X, 'COEFF VECTORS') 

WR I TE ( 6, 660 ) ( I C ( J ) , J= 1 , NC ) 

660 FORMAT !/10X. 'BASI S B ',3713) 

WRITE(6,67C> 

670 FORMAT ( ' ' ) 

DO ICO 1=1, M 
DO 90 J= 1 » NC 

XXHOLD=B( I .LEG! 1, IC( J) ) )+B! I ,LEG(2,IC( J) ) ) 
IAX(J)=XXHOLD+SlGN( .5.XXHOLD) 

90 CONTINUE 

WR I T E ( 6, 68 0 ) IB(I), KB(I), ( I AX ( J ) , J = 1 , NC ) 
680 F0RMAT!1CX,I4,2X, 12, 2X, 3713) 

100 CONTINUE 

DO 110 J=1 ,NC 

no kc ( j)=c( ic ( j ) ) 

WR I TE ( 6, 690) ( KC ( J ) , J= 1 , NC ) 

690 FORMAT (/10X, 'REL COST • , 2X , 37 I 3 ) 

WR I TE ( 6 , 640 ) KZ 

RETURN 

ENC 



C 

C 



C 

C 



SUBROUTINE SOLCHK!*) 

COMMON /TSBLK/ CO ST , C STLC ( 1 600 ) , NIC ( 1 000 ) , 

1 I C ( 1600 ) ,KC(40),KB( 120 ) , I AX ( 40 ) , 

2 I ROW! 40 ) , ICOLI 40) ,LC(40), 

3 KZ,NC»J» ICJ»L1»L2»J1»CST.XZ 
COMMON /SHARE/ C(1600),B( 80, 80), H( 80), Z, 

1 LEG( 2, 1600 ,IB( 120) ,MM,M,N 

CONVERT SUBSCRIPTS OF BASIC VARIABLES EQUAL TO ONE 
BACK TO ORIGINAL DOUBLE-SUBSCRIPTED FORM. 

1=0 

DO 20 L= 1 , MM 
10 1 = 1+1 

IF(KB(I).EQ.O) GO TO 10 
IPOW(L) = LEG( 1, I B C I ) ) 

I COL ( L )=LEG( 2, IB( I ))-MM 
20 CONTINUE 

COUNT THE NUMBER OF STEPS REQUIRED 
TO GO FROM ORIGIN BACK TO ORIGIN. 

KSTE P= 1 
K = 1 

30 DO 40 L=2 , MM 

I F ( I ROW C L I .NE. ICOL(K) ) GO TO 40 
KSTEP=KSTEP+1 

I F ( I COL ( L ) .EQ.IROW(l) ) GO TO 50 
K= L 

GC TO 30 
40 CONTINUE 

50 IFIKSTEP.LT. MM) RETURN 
CST=COST 
XZ=Z+COST 
K Z = X Z 

WRITE! 6,610) 

610 FORMAT (//1CX, • A TRAVELING SALESMAN SOLUTION IS - 
DO 60 L= 1 » MM 

WR I TE ( 6, 62 0 ) I ROW ( L ) , I COL ( L ) 

620 FORMAT (21X, • X( • , I 2, ' , ' , 12, ' ) = 1') 

60 CONTINUE 

WRITE! 6, 63 0 KZ 
630 FORMAT I/28X, 'Z =',I5) 

RETURN 1 
END 



'//) 



41 



SUBROUTINE RANKC 

COMMON / T S BL K / CO ST , C STLC I 1600 ) , NLC I l 000 ) , 

1 IC( 1600),KC<40) ,KB( 120) , I AX 1 40) , 

2 I ROW I 40 ) , I COL I 40 ) » LC I 40 ) » 

3 KZ » NC » J * ICJ»Li»L2,Jl»CST. XZ 
COMMON /SHARE/ C(1600),B( 80, 80), H( 80), Z, 

1 LEG(2,1600) »IB(120) ,MM,M,N 

NCi=NC-l 
DO 20 K= 1 , NC 1 
CM I N=C I I C I K ) ) 

LM I N=K 
K 1=K+ 1 

DO 1C L=K 1 , NC 

IF ICI IC(L) ) .GE.CMIN) GO TO 10 
CM IN=C ( I C < L) ) 

LMIN=L 
10 CONTINUE 

I CHOLD= IC ( K) 

I C (K )= ICILMIN ) 

IC(LMIN) = I CHOL D 
20 CONTINUE 
RETURN 
END 



SUBROUTINE REPARTI*) 

COMMON /TSBLK/ COST , CSTLC ( 1 600 ), NLC ( 1000 ) , 

1 ICI 1600) ,KC(40) ,K3< 120 ) . I AX (40) , 

2 I RO W I 4G ) » I CO LIAO) » LC I 40 ) , 

3 KZ,NC,J,ICJ,LI,L2,J1,CST,XZ 
COMMON /SHARE/ C( 1600),B( 80, 80), HI 80), Z, 

1 LEGI2, 1600) ,IB( 120) ,MM,M,N 

I I HOLD=0 
DO 120 11 = 1, M 

C IB IS VECTOR OF BASIC-VARIABLE SUBSCRIPTS. 

C NEGATIVE VALUES INDICATE ARTIFICIALS. 

IFIIBI II ) . GT . 0 ) GO TO 120 
DO 10 J=1 » NC 
I C J= I C I J) 

L 1=L EG 1 1 , I CJ > 

L2 = L EG 1 2 , ICJ ) 

XHLD=B 1 1 I , LI ) +BI I I ,L2 ) 

IPIV=XHLD+SI GN I .5, XHLD ) 

IF! I PIV.NE.O ) GO TO 30 
10 CONTINUE 

C IF ONE ARTIFICIAL CAN NOT BE REPLACED WITH A 
C PIVOT OPERATION, EARMARK ITS ROW FOR DELETION. 

C IF THERE IS NO PIVOT ELEMENT FOR MORE THAN ONE 
C ARTIFICIAL, TERMINATE EXECUTION ON THIS PROBLEM. 
IF! I I HOLD. NE .0 ) GO TO 20 
I I HOLD= 1 1 
GO TO 120 
20 WR I TE I 6 , 60 0 ) 

600 FORMAT I//1CX, 'CAN NOT REPLACE ARTIFICIALS') 
RETURN l 

C REPLACE ARTIFICIAL WITH A PIVOT OPERATION. 

30 I E I I I ) = I CJ 
DO 5C 1=1, M 
IFII.EQ.II ) GO TO 50 
XHLD=B I I , L 1 ) +B I I , L2 ) 

IYI J=XHLD+ SIGN! .5 , XHLD) 

IF! IYI J.EQ.O) GO TO 50 
I Q= I Y I J/ IP IV 
DO 40 L=1 , M 

40 B(I*L)=B( I ,L)-IQ*B( II »L) 

50 CONTINUE 

IF( IPIV.EQ.l ) GO TO 70 
DO 6C L= 1 , M 
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60 B( II » L ) = ( - 1 . ) * B ( I I , L ) 

70 NC=NC-1 

IF(J .EQ.NC-H ) GO TO 90 

C SHIFT NONBASIC VECTORS LEFTWARD TO FILL GAP. 
DO 80 L=J,NC 
30 I C C L > = ICIL+I ) 

C RECOMPUTE RELATIVE COSTS. 

90 DO IIC J=1,NC 
ICJ=IC( J) 

LI=LEG( I»ICJ) 

L2=LEG(2, ICJ ) 

CHOLD=CSTLC( ICJ) 

DO IOC 1 = 1, M 

IFCIBCI ).LT.O) CX=10CCOOO. 

IFCIBC IJ.GT.O) CX=CSTLC( IB ( I ) ) 

100 CHOLD=CHOLD-CX*( B( I , L 1 ) +B ( I , L2 ) ) 

110 C( I C J ) =CHO LD 
CALL RANKC 
120 CONTI NUF 

C FIND A REDUNDANT TABLEAU ROW TO DROP. 

Ml =M 
M=M— 1 

IF ( I IHOLD.NE .0 ) GO TO 150 
DO 1 AC I 1 = 1, Ml 
DO 13C J= 1 ,NC 
ICJ= IC(J ) 

L1=LEG ( 1 , ICJ ) 

L2=LEG(2, ICJ) 

XHLD=B ( II, LI )+B( II, L2) 

IYI J=XHLD+SIGN( . 5 , XHLD ) 

IF(IYIJ.NE.O) GO TO 1A0 
130 CONTINUE 
1 1 HOLD = I I 
GO TO 1 50 
1 AO CONTINUE 
RETURN 

150 I F ( I I HOLD . EQ. M 1 ) RETURN 
C SHIFT TABLEAU ROWS UPWARD TO FILL GAP. 

DO 170 I = I IHOLD, M 
I B ( I )= IB ( I +1 ) 

H ( I ) =H ( I + 1 ) 

DO 160 J= 1 , M 1 
160 B ( I, J) = B( I+1,J ) 

170 CONTINUE 
RETURN 
END 



SUBROUTINE LINCOM(*,*) 

COMMON /TSBLK/ CO ST, CSTLC ( 1600 ) ,NLC ( 1000 ) , 

1 IC ( 1600 ) ,KC( AO ) ,KB( 120 ) , I AX (AO) , 

2 IROW(AO) ,ICOL( AO) ,LC( AO) , 

3 KZ,NC, J, ICJ,L1,L2, J1,CST,XZ 

COMMON /SHARE/ C(160C),B( 80, 80). H( 80), Z, 

1 LEG( 2, 1600) ,IB( 120) ,MM,M,N 

IF ( J.LE.32 ) GO TO 20 
WRITE( 6,600) 

600 FORM AT (15X, ’CAN NOT PROCEED BECAUSE LN COM COUNT ', 
1 'WILL EXCEED 2**31 - 1 ' ) 

10 CALL TIME I TC - 1 ,TIME) 

WPITE(6»61C) 

610 FORMAT (//l 5X ,' PROGRAM LIMITATIONS FORCE ', 

1 'ALGORITHM HALT' ) 

WR ITE ( 6, 62C ) TIME 

620 FORMAT (/15X, 'ALGORITHM TIME =',F10.A,' SECONDS') 
RFTURN 2 

C INITIALIZE COUNTERS. 

20 NGEN=0 
NCST F=0 
NC NF=C 
KOST=C ST 

C INITIALIZE COMBINATION GENERATER. 
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DO 3C 1=1, J1 
30 LC( I ) = 0 

C BEGIN GENERATING. 

C J1 IS EQUAL TO J-l AND IS THE TOTAL NUMBER OF 
C NONBASIC VARIABLES TO BE TRIED IN COMBINATION 
C WITH THE J TH NONBASIC VARIABLE. 

C KT IS THE NUMBER OF BINARY NUMBERS WITH EXACTLY 
C I J 1 DIGITS WHERE THE LEFTMOST DIGIT IS ALWAYS A ONE. 

DO 170 IJ1=1,J1 
KT=2**< I Jl-1 ) 

DO 160 11=1, KT 
DO 40 1 = 1, IJ1 
I F ( LC( I ) .EQ.O) GO TO 50 
40 LC ( I ) = 0 
50 LC ( I )=1 

C INCREMENT THE NUMBER OF COMBINATIONS GENERATED BY ONE. 
NGEN=NGE N+ 1 

C COMPUTE LATEST COMBINATION'S RELATIVE COST. 

COST=C ( I CJ ) 

DO 60 1 = 1 , I J 1 

I F ( LC( II.EQ.l) COST=COST+C ( IC ( I > > 

60 CONTINUE 

IF (COST.LE .CST ) GO TO 110 
C THE FOLLOWING SEQUENCE OF STATEMENTS DOWN TO 
C STATEMENT 1 1C ARE VARIOUS CHECKS TO SEE IF 
C ENTIRE BLOCKS OF COMBINATIONS CAN BE SKIPPED. 

IF( I I.GT.l ) GO TO 70 
IF (NCNF.GT .0 ) GO TO 180 
65 WR I T E ( 6, 64 C) NGEN , KOST , NCSTF , NCNF 
640 FORMAT ( 1 5X , ' NR LN COM NECESSARY TO GENERATE =',110, 

1 / 1 5X , ' NR WITH REL COST .LE.',I4,' =',6X,I10, 

2 /I 5X , ' NR FOUND CONSTRAINT FEASIBLE =',3X,I10) 
RETURN 

70 DO 80 1 = 2, I J 1 

IF(II.EQ.2**( I-2I+1) GO TO 90 
80 CONTINUE 
GO TO 16C 
90 DO 100 1 = 1, 1 J 1 

ioo Lcm=i 

GO TO 17C 

C INCREMENT BY ONE THE NUMBER OF COMBINATIONS 
C FOUND LESS THAN OR EQUAL TO THE UPPER BOUND. 

110 NC STF = NC ST F + 1 

C CHECK THE LATEST COMBINATION FOR CONSTRAINT FEASIBILITY. 
DO 130 K=1 , M 

XHOLD=H( K) -B ( K , L 1 ) -B ( K , L2 ) 

KB ( K ) = XHOL D+S I GN ( .5,XH0LD) 

DC 1 20 I = 1 , I J 1 
IF(LCm.EQ.O) GO TO 120 

ici=icm 

XHCLD=KB(K)-B(K»LEG(1» IC I ) )-B ( K , L EG ( 2 , IC I) ) 

KB ( K ) = XHOLD+ SI GN( .5,XH0LD) 

120 CONTINUE 

IF(KB(K) .NE.0.AND.KB(K).NE.1 > GO TO 160 
130 CONTINUE 

C INCREMENT BY ONE THE NUMBER OF COMBINATIONS 
C FOUND CONSTRAINT FEASIBLE. 

NCNF=NCNF+ 1 

IF (NCNF. LE. 1000) GO TO 140 
WRITE! 6,640) NGEN , KOST , NCSTF , NCNF 
WRITE(6,650) 

650 FORMAT (15X, 'CAN NOT PROCEEO BECAUSE STORAGE IS ', 

1 'ALLOCATED FOR NO MORE THAN 1000 ', 

2 'CN FEAS LN COM' ) 

GC TO 10 

C STORE THE RELATIVE COST AND IDENTIFYING NUMBER 
C FOR THE LATEST CONSTRAINT FEASIBLE COMBINATION. 

140 CSTLC( NCNF )=COST 

NLC(NCNF)=2**( I Jl-1)+I 1-1 
I F (NCNF. EQ .1 ) GO TO 160 

C ARRANGE VALUES IN CSTLC AND NLC ACCORDING TO 
C ASCENDING ORDER OF RELATIVE COSTS. 
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NCNF1= NCNF-1 

00 150 I L= 1 , NC NF 1 
L 1 =NCNF 1- I L* 1 

IF(CSTLC(LI+1) .GE.CSTLC(LI)) GO TO 160 
CHOLO=CSTLC( LI ) 

NHOLD= NL C ( L I ) 

C STLC ( LI ) =CSTLC ( L I +1 ) 

NLC( LI )=NLC( LI +1 1 
CSTLC ( LI +1 ) = C HOLD 
150 NLC( LI +1 ) =NHOLD 
160 CONTINUE 
170 CONTINUE 

IF(NCNF.EQ.D) GO TO 65 

C CHECK EACH COMBINATION STORED ANO RANKED ABOVE 
C FOR A TRAVELING-SALESMAN SOLUTION. 

180 WRITE! 6, 64 C) NGEN , KOST , NCSTF , NCNF 
DC 2 1C 1 = 1, NCNF 

C RECOMPUTE OPTIMAL ASSIGNMENT BASIC-VARIABLE VALUES. 

DO 185 K=1,M 

XH0LD=H(K)-B(K,L1 >-B(K,L2) 

185 KB ( K ) =XHOL D+S I GN( .5»XH0LD) 

L = 0 

KM=M+1 
NL=NLC ( I ) 

C REGENERATE A COMBINATION STORED ABOVE. 

190 L = L+ 1 

LCL=MOD( NL ,2 ) 

IF(LCL.EQ.Q) GO TO 2 00 
ICL=IC(L) 

C RECOMPUTE NEW BAS IC-VAR I ABLE VALUES 
C FOR LATEST REGENERATED COMBINATION. 

DO 195 K= 1 , M 

XH0LD=KB(K)-B(K,LEG(1, ICL) )-3(K,LEG(2, ICL ) ) 

195 K B ( K ) = XHOL D+ S I GN ( .5,XH0LD) 

C AUGMENT ASSIGNMENT BASIS WITH THE NONBASIC 
C VARIABLES IN THE LATEST COMBINATION. 

KM=K M+ 1 
I B ( KM) = I CL 
KB ( KM ) = 1 
200 NL=NL/2 

1 F ( N L . GT . 0 ) GO TO 190 

C CHECK LATEST COMBINATION FOR A TRAVELING-SALESMAN TOUR. 
COST = CST LC ( I ) 

CALL S OLCHK ( £2 20) 

C AS SOON AS A TOUR IS FOUND GO TO STATEMENT 220, 

C OTHERWISE CONTINUE. 

210 CONTINUE 
RETURN 

220 RETURN 1 
END 



SUBROUTINE ASSIGN!*) 

COMMON / ASGBLK/ 0(1600), DB! 80), DUAL! 80 ) , AX ( 80), 
1 E,AD,W,INDEX,K1,JS,IR,ID,ITR 

COMMON /SHARE/ C(1600),B( 80, 80), H( 80),/, 

1 LEG(2» 1600) , IB! 120 I , MM , M, N 

WRITE! 6,35) 

35 FORMAT ( * 1 • ) 

K 1 =2 
I NDE X= 0 
Z = 0. 

NN=MM 
DO 4 1=1, M 
4 H ( I ) = 1 . 

DO 1 1=1, MM 
U=1C000. 

V= 10000. 

KA=NN* I-NN+1 
KB=NN* I 
DO 6 J=KA , KB 
IF(C(J).LT.U) U=C(J) 
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6 LEG( 1, J )=I 

IFIU.EQ.O. ) GO TO 14 
DO 27 J=KA »KB 

27 C(J)=C(J)-U 
Z= Z+U*H ( I ) 

14 K = I •*- MM 

DO 3 J=I,N,NN 
IF(CIJ).LT.V) V=C(J) 

3 LEG(2,J)=K 

IF (V .EQ.C. ) GO TO 1 
DO 28 J=I,N,NN 

28 C(J)=C(J)-V 
Z= Z+V*H{ K ) 

1 CONTINUE 

CALL L IP ( £ 33 ) 

RETURN 
33 RETURN 1 
END 



SUBROUTINE LIPI*) 

COMMON / ASGBLK / D(160P),DB( 30), DUAL! 80 ) , AX ( 80), 
1 E , AD, W, INDEX, K1,JS,IR, ID, ITR 

CGMMON /SHARE/ C(1600),B( 80, 80), H( 80), Z, 

1 LEGI2, 1600) , I B( 120) , MM,M, N 

E= .5 

DO 5 J = 1 , N 

5 D( J) =-2. 

DO 3 J = 1,M 
DC 2 I =1 , M 

2 B ( 1 1 J ) =0 • 

3 B(J,J)=1. 

DO 4 I =1 , M 
DB ( I )=-l. 

DUAL! I )=C • 

4 I B ( I ) =- I 
A D= 1 . 

I D= 1 
W=0 . 

DO 6 1=1, M 

6 W= W+H ( I ) 

I TR=0 

7 IF(ABSIW)-E) 100,100,8 

8 DO 9 J=1,N 

I F ( D ( J ) ) 10,9,9 

9 CONTINUE 
GO TO 101 

19 CX=1CC000. 

DO 11 J=1,N 
JT = J 

I F ( D ( J ) ) 12,11,11 

12 I F ( C ( J )- »0 00 1 • LT • 0 • ) GO TO 41 
IF(CX+C( J)/D(J ))11,11,13 

13 CX=-C( J)/D(J) 

JS=J 

11 CONTINUE 

DO 14 J= 1 , N 

14 C(J)=C(J)+CX*D(J) 

C( JS )=0. 

Z=Z+CX*W 

JT=J S 

DC 40 1=1, M 

40 DUAL ( I )=DUAL ( I )-CX*DB( I ) 

15 IF(W-E)1CC,100,41 

41 DX=0 • 

DO 16 J=JT,N 
I F ( C ( J )-C. 0001 ) 1 7, 17, 16 
17 C(J)=0. 

IF <D( J ) .GE .DX ) GO TO 16 
DX = D( J ) 

19 JS= J 

IF(D( JS) . EQ.-2 . ) GO TO 1 
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CONT INUE 

WRITE<6.213) JS,D(JS) 

FORMAT (lH « 3HJ S= , I4,3X»2HD=,F8.4) 

I F ( D ( J S ) ) 20 » 7 » 7 

CALL ELEM 

I F ( 10-1)528,26.26 

CALL PIVOT (C33) 

GO TO 15 
CONT INUF 
INFEA =C 
GO TO 42 5 
WRITE (6,205) 

FORM AT ( 1 H ,10HINFEASIBLE) 

INFEA =1 
RETURN 1 
CONTINUE 
1 = 1 

RFTURN 

END 



SUBROUTINE ELEM 

COMMON / ASGBLK/ 0(1600), DB( 80),DUAL( 80 ) , AX ( 80), 

E , AO , W , INDEX,K1 , JS,IR,IO, ITR 
COMMON /SHARE/ C(1600),B( 80, 80), H( 80), Z, 

LEG (2, 1600) ,IB( 120) ,MM,M,N 

I R=0 
J=JS 

DC 2 1 I = 1 , M 
AX ( I ) = 0. 

00 21 K= 1 , 2 
K2=LEG(K,J) 

AX ( I )= AX ( I ) + B ( I , K 2 ) 

21 CONTINUE 
EE=E 

DO 25 1=1, M 
IF(AX(I)-E) 25,25,810 
310 I F ( H ( I ) — E ) 82,82,23 
82 IF=I 

GO TO 13 

23 IF(IR.EQ.O) GO TO 24 
IF(ABS(H(I )- AX ( I)*X)-EE) 70,70,2 

2 I F ( H ( I )-AX(I )*X) 24,70,25 

24 X=H ( I ) /AX( I ) 

3 I R= I 

E E=E /A X ( IR ) 

GO TO 25 

70 I F ( I B ( I ) ) 74,25, 73 

73 IFUB(IR)) 25,25,4 

74 I F ( I B ( I R ) ) 4,25,3 

4 I F ( AX ( IR ) - AX ( I ) ) 25,25,3 

25 CONTINUE 
13 RETURN 

END 



SUBROUTINE PIVOT(*) 

COMMGN /ASGBLK/ 0(1600), DB( 30 ) , DUAL ( 80),AX( 80), 
1 E,AD»W, INDEX, K1,JS,IR, ID, ITR 

COMMON /SHARE/ C(16C0),B( 80, 80), H( 80), Z, 

1 LEG( 2, 1600 ) , IB( 120 ) ,MM,M,N 

84 WRITE (6,214) IR,AX(IR) 

214 FORMAT (1H+,25X,3HIR=, I 3 , 3X , 3HA X= , F 8 . 4 ) 

I B ( I R ) = J S 

AD=A D* AB S ( AX ( IR) ) 

D I = I D 

ID = INT(DI*ABS(AX( IR) ) +0 . 5 ) 

D I = I D 

IF (AX( IR ) .LT .1 . + E.AND. AX( IR) .GT. l.-E ) GO TO 1 
HX=1 . / AX ( I R) 

H ( IR )=H( IR )*HX 
DO 11 J= 1 , M 



16 

1 

213 

20 

26 

130 



101 

205 

33 

425 

523 
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11 B( IR, J )=B( IR, J)*HX 

1 DO 3 1=1, M 

IFU .EQ. IR ) GO TO 3 
I F ( ABS ( AX ( I) ).LT.E) GO TO 3 

IF(AX( I ) .LT. l.+E. AND.AX( I ) .GT.l.-E ) GO TO 4 
I F ( AX ( I) . GT . -1 .-E . AND. AX ( I ).LT .-l.+E) GO TO 5 
DO 2° J=1 , M 

29 B( I,J)=B( I f J )-AX(I)*B( IR,J) 

H ( I I =H ( I ) - AX ( I )*H( IR l 

GO TO 3 

A DO 6 J=1 , M 

6 8( I » J ) =B ( I »J)-B( I R » J ) 

H ( I ) =H ( I ) - H( IR ) 

GC TO 3 

5 DO 7 J=1,M 

7 B( I, J)=B( I ,J )+B( IR, J) 

H ( I ) = H ( I ) +H( IR ) 

3 CONTINUE 

W=W+D( JS )*H( IR ) 

DO 8 J=1 , M 

8 DB(J)=DB(J)-D(JS)*B(IR»J) 

E=0. 5/DI 

DO 90 K= 25 » 60 0 » 25 
I F ( K— I TR ) 9C , 9 1 ,95 

90 CONTINUE 
GO TO 95 

91 DO 92 J= 1 » M 

DB(J )=SIGN(AINT(DI*ABS(DB( J) )+0.5)/DI ,DB( J) ) 

DO 92 1=1, M 

92 B( I, J)=SIGN(AINT(DI*ABS(B( I , JI)+0.5)/DI,B( I, J)l 
DO 93 1=1, M 

93 H ( I ) =S IGN( AINT(DI*ABS(H( I )) +0 . 5 ) / D I , H (II ) 

W=AINT (DI*W+0. 5I/DI 

95 DO 34 J= 1 , N 

D ( J I =0 . 

DO 33 K=1,K1 
K2 = L EG ( K , J ) 

I F ( K2 . EQ.CI GO TO 33 
D(JI=D(J)+DB(K2) 

33 CONTINUE 
IF(ABS(D(J))-EI35,35,34 

35 D ( J I =0 . 

34 CONTINUE 

I TP= I TR+1 

WRITE (6,210) ITR,W,Z, ID, AD 
210 FORMAT ( 1H , 5X ,4H I TR= , I 5 , 5X , 2HW= , F 1 2 . 8 , 5X , 2HZ= , 

1 F12.5, 5X,3HID=,I6,3X,3HAD= ,F12.4I 

I F ( ITR.EQ.300IGO TO 21 
RETURN 

21 WRITE( 6,220 

220 FORMAT (//6X, ' ITERATION LIMIT EXCEEDED') 

RETURN 1 
END 



SUBROUTINE R AND ( * ) 

COMMON /SHARE/ C( 1600),B( 8>0, 80), H( 80), Z, 
1 LEG(2, 1600) , IB( 120) ,MM,M,N 

DATA 1/3125/ 

DO 10 J=1,N 
I=MOD(3125*I ,2048) 

X=( 100. *1 ) /2048. + .5 
K= X 

C( J)=K 
10 CONTI NUF 
MM1=MM + 1 
DO 20 J= 1 » N, MM1 
C( J)=999. 

20 CONTINUE 
RETURN 1 
END 
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SUBROUTINE T I ME I T ( N , T I HE ) 



c 

r 


N=< 


D STARTS 


CLOCK; N=- 1 STOPS CLOCK 




IT 


=N +2 






GC 


TO ( 20 , 


10 ) , IT 


10 


CALL TIMON(M) 




T I MEM=M 






RETURN 




20 


CALL TIMOFF(M) 




T I ME=M 






T IME=(TIMEM- TIME)*. 000026 




RETURN 






END 




C PLACE 


A /* CARD HERE. 


//ASM. 


SYSIN DD * 




TIMEALL 


CSECT 








ENTRY 


TIMON, TIMOFF 


TIMON 




SAVE 


(14,12) ENTRY VIA 






USING 


TIMON, 12 






LR 


12,15 






ST 


13, TEMPI 






LA 


13, SAVE 1 






L 


2 , 0 ( 1 , 0 ) 






L 


3, TOTIME 






ST 


3, CLOCKR 






ST 


3, 0(2,0) 






STIMER 


TASK,TUINTVL=CLOCKR 






L 


13, TEMPI 



TIMON(N)- 15.5 



EXIT 

TIMOFF 



TOTIME 
CLOCKR 
SAVE 1 
TEMPI 



RETURN 

SAVE 

USING 

LR 

ST 

LA 

L 

TTIMER 

ST 

L 

RETURN 
CNOP 
DC 
DS 
DS 
DS 
END 
A /* 



114,12) ,T,RC=0 

(14.12) ENTER 
TIMOFF, 12 

12,15 
13, TEMPI 
13 , S AV E 1 
2 » 0 < 1 , 0 ) 

CANCEL 

C, 0(2,0) 

13, TEMPI 

(14.12) ,T,RC=0 
0,4 

X* 7FFFFFFF • 

F 

18F 

F 



VIA -CALL TIMOFF(N)- 



C PLACE 
//GO . FTC 6F00 1 
/ /GO . SYS I N DD 
C PLACE DATA 
C PLACE AN 



HERE. 

S PAC E = ( C YL , ( 10,1) > 



CARD 
DD 

* 

CARDS HERE. 
ORANGE /* CARD 



HERE, 
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APPENDIX B COMPUTER OUTPUT FOR A SELECTED PROBLEM 



2 

COST MATRIX 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


999 


24 


18 


22 


31 


19 


33 


25 


30 


26 


2 


15 


999 


19 


27 


26 


32 


25 


31 


28 


18 


3 


22 


23 


999 


23 


16 


29 


27 


18 


16 


27 


4 


24 


31 


18 


999 


19 


13 


28 


9 


19 


27 


5 


23 


18 


34 


20 


999 


31 


24 


15 


25 


8 


6 


24 


12 


17 


15 


10 


999 


11 


16 


21 


31 


7 


28 


15 


27 


35 


19 


18 


999 


21 


21 


19 


8 


13 


24 


18 


13 


13 


22 


25 


999 


29 


24 


9 


17 


21 


18 


24 


27 


24 


34 


31 


999 


18 


10 


18 


19 


29 


16 


23 


17 


18 


31 


23 


999 



OPTIMAL ASSIGNMENT SOLUTION 

BASIS B 

82 0 
11 1 
75 1 

38 1 

50 1 

29 1 

3 0 

94 1 

6 1 
81 0 
74 0 

62 1 
83 1 

4 0 

55 0 

69 0 

57 1 

48 0 

90 0 

z = 140 



Alternate solution procedures for this ten-city problem are 
discussed in References 8 and 9, 
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TRY X( 8) = 1 REL COST = 0 



TRY X(97) = 1 REL COST = 1 
TRY THE (2** 1) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 1 

NR WITH REL COST .LE. 2 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(68) = 1 REL COST = 2 
TRY THE (2** 2) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 2 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(99) = 1 REL COST = 2 
TRY THE (2** 3) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 2 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(20) = 1 REL COST = 2 
TRY THE (2** 4) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 2 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(84) = 1 REL COST = 2 
TRY THE (2** 5) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 3 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(20) AND 

TRY THE (2** 4) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 3 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 3 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 3 

NR WITH REL COST .LE. 3 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 
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RETURN TO X(97) AND 

TRY THE (2** 1) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 1 

NR WITH REL COST , LE. 3 = 1 

NR FOUND CONSTRAINT REASIBLE = 0 

TRY X(58) = 1 REL COST = 3 
TRY THE (2** 6) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 3 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X( 2) = 1 REL COST = 3 
TRY THE (2** 7) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 3 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X( 9) = 1 REL COST = 3 
TRY THE (2** 8) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 3 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(19) = 1 REL COST = 3 
TRY THE (2** 9) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 3 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(13) = 1 REL COST = 3 
TRY THE (2**10) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 3 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(52) = 1 REL COST = 3 
TRY THE (2**11) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 3 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(65) = 1 REL COST = 3 
TRY THE (2**12) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(52) AND 

TRY THE (2**11) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 
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RETURN TO X(13) AND 

TRY THE (2**10) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 1 

RETURN TO X(19) AND 

TRY THE (2** 9) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X( 9) AND 

TRY THE (2** 8) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X( 2) AND 

TRY THE (2** 7) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(58) AND 

TRY THE (2** 6) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(84) AND 

TRY THE (2** 5) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 12 
NR WITH REL COST .LE. 4 = 9 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(20) AND 

TRY THE (2** 4) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 9 

NR WITH REL COST .LE. 4 = 7 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 6 

NR WITH REL COST .LE. 4 = 5 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 3 

NR WITH REL COST .LE. 4 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 
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RETURN TO X(97) AND 

TRY THE (2** 1) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 1 

NR WITH REL COST .LE. 4 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(28) = 1 REL COST = 4 
TRY THE (2**13) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 4 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(96) = 1 REL COST = 4 
TRY THE (2**14) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 4 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(17) = 1 REL COST = 4 
TRY THE (2**15) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 4 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(92) = 1 REL COST = 4 
TRY THE (2**16) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 5 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(17't AND 

TRY THE (2**15) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 5 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(96) AND 

TRY THE (2**14) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 5 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(28) AND 

TRY THE (2**13) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 5 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(65) AND 

TRY THE (2**12) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 
NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 
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RETURN TO X(52) AND 

TRY THE (2**n) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 

NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(13) AND 

TRY THE (2**10) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 

NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 1 

RETURN TO X(19) AND 

TRY THE (2** 9) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 

NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X( 9) AND 

TRY THE (2** 8) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE =16 
NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X( 2) AND 

TRY THE (2** 7) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 

NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(58) AND 

TRY THE (2** 6) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 15 

NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(84) AND 

TRY THE (2** 5) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 17 

NR WITH REL COST .LE. 5 = 15 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(20) AND 

TRY THE (2** 4) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 12 

NR WITH REL COST .LE. 5 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 7 

NR WITH REL COST .LE. 5 = 7 

NR FOUND CONSTRAINT FEASIBLE = 0 
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RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 3 

NR WITH REL COST .LE, 5 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(97) AND 

TRY THE (2** 1) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 1 

NR WITH REL COST .LE. 5 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(86) = 1 REL COST = 5 

TRY THE (2**17) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 5 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(54) = 1 REL COST = 5 
TRY THE (2**18) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 5 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(66) = 1 REL COST = 5 
TRY THE (2**19) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 5 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(71) = 1 REL COST = 5 
TRY THE (2**20), - 1 Lt( COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 5 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(25) = 1 REL COST = 5 
TRY THE (2**21) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 2 

NR WITH REL COST .LE. 5 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 

TRY X(85) = 1 REL COST = 5 
TRY THE (2**22). - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 6 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(25) AND 

TRY THE (2**21) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 6 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 
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RETURN TO X(71) AND 

TRY THE (2**20) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 6 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(66) AND 

TRY THE (2**19) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 6 = 3 

NR FOUND CONSTRAINT FEASIBLE = 1 

RETURN TO X(54) AND 

TRY THE (2**18) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 6 = 3 

NR FOUND CONSTRAINT FEASIBLE = 1 

RETURN TO X(86) AND 

TRY THE (2**17) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 4 

NR WITH REL COST .LE. 6 = 3 

NR FOUND CONSTRAINT FEASIBLE = 1 

RETURN TO X(92) AND 

TRY THE (2**16) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 
NR WITH REL COST .LE. 6 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(17) AND 

TRY THE (2**15) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE =16 
NR WITH REL COST .LE. 6 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(96) AND 

TRY THE (2**14) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 
NR WITH REL COST .LE. 6 = 11 

NR FOUND CONSTRAINT FEASIBLE = 1 

A TRAVELING SALESMAN SOLUTION IS — 

X( 2, 1) = 1 
X( 8, 5) = 1 
X( 4, 8) = 1 
X( 5,10) = 1 
X( 3, 9) = 1 
X( 1, 3) = 1 
X( 7, 2) = 1 
X( 6, 7) = 1 
X(10, 6) = 1 
X( 9, 4) = 1 

Z = 146 
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RETURN TO X(28) AND 

TRY THE (2**13) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 16 

NR WITH REL COST .LE. 6 = 11 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(65) AND 

TRY THE (2**12) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 40 

NR WITH REL COST .LE. 6 = 31 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(52) AND 

TRY THE (2**11), - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 37 

NR WITH REL COST .LE. 6 = 29 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(13) AND 

TRY THE (2**10) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 34 

NR WITH REL COST .LE. 6 = 27 

NR FOUND CONSTRAINT FEASIBLE = 1 

RETURN TO X(19) AND 

TRY THE (2** 9) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 31 

NR WITH REL COST .LE. 6 = 25 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X( 9) AND 

TRY THE (2** 8) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 28 

NR WITH REL COST .LE. 6 = 23 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X( 2) AND 

TRY THE (2** 7) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 25 

NR WITH REL COST .LE. 6 = 21 

NR FOUND 'CONSTRAINT FEASIBLE = 0 

RETURN TO X(58) AND 

TRY THE (2** 6) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 22 

NR WITH REL COST .LE. 6 = 19 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(84) AND 

TRY THE (2** 5) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 31 

NR WITH REL COST .LE. 6 = 21 

NR FOUND CONSTRAINT FEASIBLE = 0 
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RETURN TO X(20) AND 

TRY THE (2** 4) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 15 

NR WITH REL COST .LE. 6 = 13 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(99) AND 

TRY THE (2** 3) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 7 

NR WITH REL COST .LE. 6 = 7 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(68) AND 

TRY THE (2** 2) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 3 

NR WITH REL COST .LE. 6 = 3 

NR FOUND CONSTRAINT FEASIBLE = 0 

RETURN TO X(97) AND 

TRY THE (2** 1) - 1 LN COM 

NR LN COM NECESSARY TO GENERATE = 1 

NR WITH REL COST. LE. 6 = 1 

NR FOUND CONSTRAINT FEASIBLE = 0 



**Z = 146 IS OPTIMAL ** 

ALGORITHM TIME = 1.5043 SECONDS 
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